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Abstract. We consider the n — » limit of the DGLAP splitting function "f gg (a) at all orders in 
the strong coupling constant, a, by analysing the leading order large Nf form of the associated 
c?-dimensional critical exponent. We show that for unpolarized scattering the pole at n = 
which appears in successive orders in perturbation theory is absent in the resummed expression. 
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The DGLAP equation, jy], is widely used to evolve the parton structure functions of the 
nucleon constituents over a large range of energy scales. Central to the equation are the splitting 
functions which depend on the variable x which represents the fraction of the momentum carried 
by that parton in the nucleon. Alternatively one can work with the anomalous dimension of the 
underlying twist- 2 operators built out of the quark and gluon fields of QCD which depend on the 
variable n. They are related to the splitting functions via a Mellin transform restricted to the unit 
interval in x. It is generally accepted that the solution of the DGLAP evolution which effectively 
represents perturbative QCD, fits the data extremely well. See, for example, @. Moreover, it 
appears to transcend the region where the perturbative approximation ought not to be valid. 
Therefore, one would hope that either by resumming perturbation theory in some fashion or 
developing non-perturbative methods, such as the BFKL formalism Q, it might be possible to 
begin to explore phenomena in the more extreme non-perturbative regions. One recent issue, 0], 
has been the problem of understanding the n — > behaviour of the DGLAP splitting functions 
in QCD. It is known that in perturbation theory each term of the expansion of the operator 
anomalous dimensions in the strong coupling constant has successively higher order poles at 
n = 0, H 2> 01- However, other considerations suggest that the function is finite at n = 0, @, §], 
though there is some disagreement with this point of view, ||. In this letter we provide some 
more insight into the n — > limit of the gluon-gluon splitting function, r y gg (a), where a = a s / (47r) 
is the strong coupling constant. This is achieved by examining the large Nf result for 7 95 (a) 
which has been computed in [jl^] where Nf is the number of quark flavours. Essentially the 1/Nf 
expansion sums a different set of Feynman diagrams from those which are ordinarily computed 
in the loop or coupling constant expansion of perturbation theory. Moreover, the resummation 
is to all orders in a. This technique has been used in [Tl| , 12, 1C, [b|] to determine the anomalous 
dimensions of all the twist-2 operators used in deep inelastic scattering for both unpolarized and 
polarized processes at 0(1/Nf). The all orders results are expressed as a function of d, where d 
is the dimension of spacetime, known as critical exponents. The coefficients of the corresponding 
renormalization group function are deduced from knowledge of the d-dimensional fixed point of 
the QCD /3-function and properties of the critical renormalization group equation. The results 
have been shown to be in agreement with all known explicit two and three loop perturbative 
results, ]|, |(| 0, [lj], 15 1, which puts the validity of the large Nf results in relation to deep inelastic 
scattering on a firm footing. Hence, in this letter we will examine the n — > limit of the singlet 
gluon splitting function at 0(1/Nf) and all orders in a. 

First, we recall the basic formalism of the problem. As the flavour singlet gluonic twist-2 
operator mixes with the fermionic operator under renormalization in perturbation theory, one 
has to deal with a matrix of anomalous dimensions, 7«(ct). Since we will be considering it in 
the large Nf expansion we define the coefficients of its perturbative expansion formally by 




where 

l qq (a) = aia + (a 2 %Nf + a 22 )a 2 + (a 3 i-/V/ + a 32 N f + a 33 )a 3 + 0(a 4 ) 

l gq {a) = ha + (b 21 N f + b 22 )a 2 + (631JV/ + 632^/ + 6 33 )a 3 + 0(a 4 ) 

l qg {a) = ciNfa + c 2 Nfa 2 + (c 31 N f 2 + c 32 N f + c 33 )a 3 + 0(a 4 ) 

7g9 (a) = (d u N f + d 12 )a + (d 21 N f + d 22 )a 2 + {d 31 N f 2 + d 32 N f + d 33 )a 3 + 0(a 4 ) . (2) 

The explicit values of the coefficients to two loops as a function of n are given in || ^, [?J and 
exact values for the low operator moments at three loops are found in 14, 15 1. In particular we 
note 

dn = \t{R) (3) 
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where T(R) is given by tr(T a T b ) 
structure constants are f abc . In 



T(R)5 ab and T a are the generators of the colour group whose 

[| it was argued that the set of coefficients of the C2 (G) sector 

C 2 {G) 



of 7 99 (a) at leading order in 1/iVf, corresponding to d^ at the Ith. loop, could be written 
compactly in the critical exponent form as 



vC 2 (G) 



(«c) 



[32fi b n z + 32^n + 32// + 8//n 4 + 16//n 3 



120/tV 



128/i 4 n 



160/i 4 



32/t 3 n 4 



64/x'V + 160/i 3 n 2 + 192/Tn + 316/i 3 



306/i 2 



+ 48/Tn 4 + 96/Tn" 3 - 78/Tn z - 126/i z n 
+ 31/xn + 146/i + 7n 4 + 14n 3 + 7n 
/[8n(n - l) 3 (n + 2)(n 2 - l)r(2 - (i)T(jJi + n)\ 



4 , r ^3 , v ^2 28]r(n + 2 



31/m 4 - 62/in 3 

- M)r(M) 



[32/u 5 n 2 + 32/i 5 n + 32// 



144/xV 



144//n - 160/i 4 



3„4 



8//n 3 + 240/i 3 n 2 + 244// 3 n + 316/i 3 + 16/i 2 n 4 + 32/i 2 n 

,n 3 

6n 2 - 14n - 28] 



4/Tn 
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180/i 2 n 2 - 196// 2 n - 306/i 2 - 20/m 4 - 40/m 3 + 59/m 2 



+ 79/xra + 146^ + 8n 4 + 16n 3 
/[8n(/i- l) 3 (n + 2)(n 2 - 1)] 



+ 2(/i-l)5i(n) 



(2ii-l)(jM-2)T(R) 



where 5j(n) = E~ 1 l/r l ,d = 2/i, / acrf / M = C 2 (G)<5 a6 and 



n? 



(2/x - l)Qx - 2)r(2/x) 

4r2(^)r(^ + i)r(2- M ) 



(4) 



(5) 



We recall that a feature of the large Nf approach to computing information on the perturbative 
coefficients of (||) was that the anomalous dimensions of the eigen-operators of ([!]) at criticality 
were determined and denoted by A±(a c ) = ^±,i( a c)/N} , flTC|]. That eigen-operator which 
was predominantly gluonic in content corresponds to the eigen-critical exponent A+(a c ) whilst 
A_(a c ) corresponds to the dimension of the mainly fermionic eigen-operator. The location of 
the fixed point, a c , is given by the non-trivial zero of the d-dimensional QCD /3-function and in 
the present notation, @, [L7|, |l|, 



3e 



AT(R)Nf 
+ 



+ 



T 2 (R)Nj 



™C A e 
16 



— Cp H Ca 

16 16 



64 128 



64 128 1 



1024 



[(288C(3) + 214)C7 F + (480C(3) - 229)C A ] e 5 + 0(e 6 



to 0(1 /Nf) in the large Nf expansion where d 



From (0) and (||), we find, for instance, 

128[9n 4 + 18n 3 + 79n 2 + 70n + 32]5 3 (n) 



(6) 



4 — 2e and £(r) is the Riemann zeta function. 
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243(n + 2)(n+ l) 2 (n- l)n 2 
256 [9n 4 + 18n 3 + 79n 2 + 70n + 32]S 3 (n) 



243(n + 2)(n + l) 2 (n-l)ra 2 
128[63n 6 + 189n 5 + 821n 4 + 1327n 3 + 1176n 2 + 544n + 96]S?(n) 
243(n + 2)(n + l) 3 (n- l)n 3 
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- 256 [4n 10 + 20n 9 + 17n 6 



52n' 



414n 6 - 976n 5 - 1521n 4 - 1516n 2 



- 930n 2 - 320n - 48] Si (n)/ [243 (n + 2)(n + l) 4 (n - 1)« 4 ] 
2048[n 2 + n + 1]C(4) j 1024Si(re)C(4) 102405i(n)C(3) 
~ 27(re + 2)(n 2 - l)n 27 243 

128[3n 6 + 9n 5 + 307n 4 + 599n 3 + 746n 2 + 448n + 96]C(3) 
+ 243(n + 2)(n + l) 2 (n-l)n 2 

+ 4[155re 12 + 930n n + 1455n 10 - 1250n 9 - 9879n 8 - 21786n 7 

- 47107n 6 - 80550re 5 - 97392n 4 - 78336n 3 - 40000n 2 

- 11776re - 1536]/[243(n + 2)(n + l) 5 (re - l)n 5 ] . 
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This illustrates earlier remarks in that this coefficient clearly has a fifth order pole at n = 
which is a degree larger than the previous loop order, [jTc( |. Clearly to all orders this degree of 
divergence will increase. However, since encodes the structure of this particular coefficient of 
lgg{ a ) to an orders in the coupling constant in leading order in 1/Nf as a function of re, we can 
examine the form of the exponent for re = 0. In other words prior to making the connection with 
explicit perturbation theory. Therefore setting n = in (|j) we can determine if the expression is 
divergent which would indicate that such a pole persists in J gg (a) or if it is finite. The latter case 
would indicate consistency with the supposition from more general principles, |4], Qj. First, by 
inspection of (Q) the Si (re) independent terms each have simple poles in re. However, by careful 
examination of the numerators of each term and in particular their re- independent terms, it is 
evident that the residue at n = is the same for both whilst the relative sign between the terms 
means that overall the residue is zero for the simple pole at n = 0. Also for the term involving 
Si (re) one writes 



Si+i{n) 



^ (/) (n + l) 



(1) 



(8) 



for general /, which vanishes in the limit n — > where ip{x) is the derivative of the logarithm of 
the Euler T-function. Hence the exponent (0) is in fact finite at n = 0. This is consistent with the 
resummed splitting function being non-singular at this point. However, it is important to recall 
that we have only demonstrated this property for that part of the series which is leading order 
in the 1/Nf expansion. Moreover that this part is not inconsistent with general considerations 
is, we believe, an important observation and the central result of this article. Furthermore, we 
record that when n = the exponent (|j) is 



A 



c 2 (G) 



iflc) 



n=0 



[16// - 64// + 94/i 2 - 59/i + 14] - - 2) [00 - 1) - V(3 - y)] 



+ fi(8fi 4 - 26// + 23// - 4) 



8(2 M -1)(m-1) 3 (m-2) s 



(9) 



For completeness, we also comment on the n — > limit of the critical exponents of the other 
twist- 2 operators at 0{\/Nf). We recall that, Q, 



A-,i(a c ) 



My-1)C 2 (R) V ° 
(2 fi -l)Q l -2)T(R)N f 



Q-l)(re-l)(2/t + n-2) 
2fi(n + n- l)(/i + re-2) 



+ h/>(/z - 1 + re) - V(m)] 



r(n-i)r(2/x) 



4(// + re - l)(u + re - 2)r(2/x - 1 + re) 



^n 2 + re + 2/i - 2) 2 + 2(// - 2)(re(re - 1)(2// - 3 + 2re) + 2(u - 1 + re)) 



(10) 
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for the flavour singlet fermionic operator where T a T a = C2(R)- The critical exponent for the 
non-singlet operator dimension corresponds to the first two terms of (|io| ) and it is clear that in 
the limit we are interested in that that exponent is finite. In particular 



non-sinslet / 

A_ (1 5 (a c 



4(2^ 2 -4/i + l)C 2 (fiK 
«=o (2// - - 2) 2 T(R)N f ' 1 ' 



By contrast, the singlet fermionic eigen-operator, which in perturbation theory corresponds to 
the combination of coefficients (an — bi\c\/d\\) at lih loop, diverges at n = having the behaviour 

A - l(ac) n( P - 2YT(R)N f (12) 

as n — > with the singularity arising from the final term of fllQ| ) which corresponds to Feynman 
diagrams with the operator inserted in a closed quark loop. 

To conclude with we have demonstrated that for unpolarized scattering the behaviour of 
7 55 (a) at n = at leading order in the large Nt expansion is consistent with the function being 
finite at this point. However, it is worth putting this result in context with other observations 
from the BFKL formalism which has also been used to study the anomalous dimensions discussed 
here. In |2lj] it was shown that the C2(G) sector of 7 ff9 (o) in the BFKL approach not only does 
not contain any poles in n but has in fact three branch points. Indeed this does not appear to 
be just a feature of QCD. If one examines the toy model of scalar </> 3 theory in six dimensions 
the complete analogous anomalous dimension was resummed at leading order without reference 
to the small x limit and again only branch points and no poles were observed, |^]. For other 
anomalous dimensions in QCD, the non-singlet and polarized dimensions also do not have poles 



at n = but branch points, p3]. (For a recent review of these issues see, for example, [24].) 
Therefore it is reassuring that our large Nf observation is also not inconsistent with the BFKL 
formalism. In light of this it would be interesting to go beyond our leading order 1/Nf analysis 
to confirm the absence of the n = pole at next order as well as being able to understand the 
branch point structure in the large Nf point of view. To compute the 0(1/ Nf) correction to 
A+(a c ), though, would involve the evaluation of a large set of 1/Nf graphs which contain, for 
example, several six loop diagrams and is therefore, we believe, not attainable in the near future. 
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